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Abstract 

We study the critical properties of the weakly disordered p-component random 
Heisenberg ferromagnet. It is shown that if the specific heat critical exponent of the 
pure system is positive, the traditional renormalization group (RG) flows at dimensions 
D = 4 — e, which are usually considered as describing the disorder-induced universal crit- 
ical behavior, are unstable with respect to replica symmetry breaking (RSB) potentials 
as found in spin glasses. It is demonstrated that the RG flows involving RSB potentials 
lead to fixed points which have the structure known as the 1 step RSB, and there exists 
a whole spectrum of such fixed points. It is argued that spontaneous RSB can occur 
due to the interactions of the fluctuating fields with the local non-perturbative degrees 
of freedom coming from the multiple local minima solutions of the mean-field equations. 
However, it is not clear whether or not RSB occurs for infinitesimally weak disorder. 
Physical consequences of these conclusions are discussed. 
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1 INTRODUCTION 



It has been many years since the effects produced by weak quenched disorder 
on the critical phenomena were considered to be qualitatively understood. In 
the most general terms the traditional point of view could be summarized as 
follows. 

If the disorder is weak (e.g. the concentration of impurities is small), its 
effect on the critical behavior in the vicinity of the phase transition point T c 
remains negligible so long as the correlation length R c is not too large, i.e. for 
temperatures T not too close to T c . In this regime the critical behavior will be 
essentially the same as in the pure system. However, as r = (T — T c )/T c — > 
and R c {t) becomes larger than the average distance between impurities, their 
influence can become crucial. 

As T c is approached the following change of length scale takes place. First, 
the correlation length of the fluctuations becomes much larger than the lattice 
spacing (which we take to be unity), and the system "forgets" about the lat- 
tice. The only relevant scale that remains in the system in this regime is the 
correlation length R c (t). Then, in the close vicinity of the critical point, R c 
grows and becomes larger than the average distance between the impurities, 
so that the effective concentration of impurities, measured with respect to the 
correlation length, becomes large. Such a situation is reached for an arbitrary 
small concentration u of impurities. The strength of disorder, as scaled by u, 
affects only the width of the temperature region near T c in which the effective 
concentration gets large. If uR® 3> 1, where D is the spatial dimensionality, 
one has no grounds, in general, for believing that the effect of impurities will 
be small. 

A very simple general criterion has been discovered, the so-called Harris cri- 
terion [Q], which makes it possible to predict the effect of impurities qualitatively 
from only the critical exponents of the pure system. According to this criterion 
the impurities change the critical behavior only if a, the specific heat exponent 
of the pure system, is greater than zero (i. e. the specific heat of the pure 
system is divergent at the critical point); more properly, this criterion should 
be stated as Dv > 2, where v is the correlation length exponent.) According 
to the traditional point of view, when this criterion is satisfied, a new universal 
critical behavior, with new critical exponents, is established sufficiently close to 
the phase transition point @, Q . These new exponents were found to satisfy the 
opposite criterion (i. e. a was negative) and were therefore apparently stable. 
In contrast, when a < (the specific heat is finite), the impurities appear to 
be irrelevant, i.e. their presence does not affect the critical behavior. 

We now consider this point in more detail. Near the phase transition point 
the D-dimensional Ising-like systems can be described in terms of the scalar 
field Ginsburg-Landau Hamiltonian with a double-well potential: 
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H= I d D x 



^(V0(x)) 2 + 1 -[r- 5r(x)]<f(x) + \g^\x) 



XI) 



Here the quenched disorder is described by random fluctuations of the effec- 
tive transition temperature 8t(x) whose probability distribution is taken to be 
symmetric and Gaussian: 

P[6t] =p exp^-j-Jd D x(5r(x)) 2) j , (1.2) 

where u <C 1 is the small parameter which describes the disorder, and p$ is the 
normalization constant. In Eq. fll.lQ t ~ (T — T c ) and for notational simplicity, 
we define the sign of 8t(x) so that positive fluctuations lead to locally ordered 
regions, whose effects are the object of our study. 

Configurations of the fields <p(x) which correspond to local minima in H 
satisfy the saddle-point equation: 

- A(/)(x) + T<p(x) + g<f{x) = 6r(x)(f)(x) . (1.3) 

Such localized solutions exist in regions of space where r — 5t(x) assumes neg- 
ative values. Clearly, the solutions of Eq. Ql.3|) depend on a particular configu- 
ration of the function 8r{x) being inhomogeneous. Let us estimate under which 
conditions the quenched fluctuations of the effective transition temperature are 
the dominant factor for the local minima field configurations. 

Let us consider a large region of a linear size L » 1. The spatially 
average value of the function St(x) in this region could be defined as follows: 

6r(n L ) = 4? / d D x8r(x) . (1.4) 



L 



Correspondingly, for the characteristic value of the temperature fluctuations 
(averaged over realizations) in this region one gets: 



5r L = ^St 2 (Q l ) = V^iL- ' 2 . (1.5) 

Then, the average value of the order parameter 4>{£Il) in this region can be 
estimated from the equation: 

-T + g<f = 5T(n L ) . (1.6) 

One can easily see that if the value of r is sufficiently small, i. e. if 

5r(tt L ) » r (1.7) 

then the solutions of Eq. (|1.6|) are defined only by the value of the random 
temperature: 

0(fi L ) ~ ±(MM)V2 . (1 . 8) 
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Now let us estimate up to which sizes of locally ordered regions this may occur. 
According to Eq. (|OD the condition 5t(Q l ) >> r yields: 



u l ' D 



L«~ m5 . (1.9) 

On the other hand, the estimation of the order parameter in terms of the 
saddle-point equation ( |1.6|) could be correct only at scales much larger than 
the correlation length R c ~ t~ v . Thus, one has the lower bound for L: 

L»t- u . (1.10) 

Therefore, quenched temperature fluctuations are relevant when 

or 

t 2 - uD « u . (1.12) 

According to the scaling relations, one has 2 — vD = a. Thus one recovers the 
Harris criterion: if the specific heat of the pure system is positive, then in the 
temperature interval, 

r<r* =u 1/a (1.13) 

the disorder becomes relevant. This argument identifies l/« as the cross-over 
exponent associated with randomness. 

The above considerations also demonstrate another very important point: 
in the disorder dominated region one finds a macroscopic number of local mini- 



mum solutions of the saddle-point equation ( |1.3| ). Indeed, for a given realization 
of the random function 5r{x) there exists a macroscopic number of spatial "is- 
lands" where r — 5r(x) is negative (so that the local effective temperature is 
below T c ), and in each of these "islands" one finds two local minimum configura- 
tions of the field: one which is "up", and another which is "down", as indicated 
in Eq.( |1.8| ). These local minimal energy configurations are separated by finite 
energy barriers, whose heights become larger as the size of the "islands" are 
increased. 

Now, if one is interested in the critical properties of the system, one has to 
integrate over all local field configurations up to the scale of the correlation 
length. This type of calculation is usually performed using a Renormalization 
Group (RG) scheme, which self-consistently takes into account all the fluctua- 
tions of the field on length scales up to R c . 

The point, however, is that the traditional RG approach is only a perturba- 
tive theory (albeit a powerful one) in which one treats the deviations of the field 
around the ground state configuration, and it can not take into account other 
local minimum configurations which are "beyond barriers". This problem does 
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not arise in the pure systems, where the solution of the saddle-point equation 
is unique. However, in a situation like that discussed above, when one gets 
numerous local minimum configurations separated by finite barriers, the direct 
application of the traditional RG scheme may be questioned. 

In a systematic approach one would like to integrate in an RG way over 
fluctuations around the local minima configurations. Furthermore, one also 
has to sum over all these local minima up to the scale of the correlation length. 
In view of the fact that the local minima configurations are defined by the 
random quenched function Sr(x) in an essentially non-local way, the possibility 
to implement successfully such a systematic approach seems rather hopeless. 

On the other hand there exists another technique which has been devel- 
oped specifically for dealing with systems which exhibit numerous local minima 
states. It is the Parisi Replica Symmetry Breaking (RSB) scheme which has 
proved to be essential in the mean-field theory of spin-glasses (see e.g. 0). 
Recent studies show that in certain cases the RSB approach can also be gen- 
eralized for situations where one has to deal with fluctuations as well 
I- 

In this paper we are going to study the critical properties of weakly disor- 
dered systems in terms of the RG approach generalized to take into account 
the possibility of the RSB phenomena. The model we will treat is the ferro- 
magnetic 0(n) or Heisenberg model with weak random interactions. The idea 
is that hopefully, like in spin-glasses, this type of generalized RG scheme self- 
consistently takes into account relevant degrees of freedom coming from the 
numerous local minima. In particular, the instability of the traditional Replica 
Symmetric (RS) fixed points with respect to RSB indicates that the multiplicity 
of the local minima can be relevant for the critical properties in the fluctuation 
region. 

It will be shown in the next Section that, whenever the disorder appears to be 
relevant for the critical behavior (in accordance with the Harris criterion), the 
usual RS fixed points (which used to be considered as providing new universal 
disorder-induced critical exponents) are unstable with respect to "turning on" 
an RSB potential. In the presence of such a potential the RG flows actually 
go to another type of fixed point having a structure known as 1-step RSB. 
The 1-step RSB structure is described by one parameter x (0 < x < 1), 
which is the coordinate of the RSB step, and this parameter remains arbitrary 
within the framework of the RG scheme. Therefore, within the framework of 
the formal RG consideration one finds a whole line of the fixed points (instead 
of the unique fixed point in the RS subspace), and correspondingly one obtains 
a whole spectrum of critical exponents. 

Formally, the value of the parameter xq at the fixed point is defined by 
the "initial conditions" for the RG equations. This situation is qualitatively 
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different from the traditional RS one, where the fixed point appears to be 
universal, and does not depend on the the details of the starting values of the 
parameters of the Hamiltonian. Here, the existence of such RSB fixed points 
indicates that the actual values of the critical exponents could be non-universal, 
being dependent on the concrete characteristics of the disorder involved. 

In Section 3 we discuss the problem of the "initial conditions" for the RG 
calculations in more detail. The crucial problem for the present approach is 
to understand whether RSB is inherent in the random bond model or not. If, 
for instance, RSB does NOT occur spontaneously for the weakly random bond 
Heisenberg model, then one remains in the replica symmetric subspace of poten- 
tials and the traditional approach^], [| is probably valid. In contrast, it might 
be that disorder is always accompanied by a small amount of RSB. In that 
case, critical exponents, amplitude ratios and the like, would be determined 
by the RSB fixed point we find here. Based on general physical arguments 
we propose a mechanism whereby the multiple local minimum solutions dis- 
cussed above can provide RSB interactions between the fluctuating fields. It 
is argued that to sum over all the local degrees of freedom, one has to sum 
over discrete local minima solutions first, and then one could initiate the RG 
integrations over the fluctuating fields. The point is that the fluctuating fields 
themselves are the deviations from the local minimum configurations, and that 
is why the summation over these quenched discrete degrees of freedom (to- 
gether with the replica averaging over the quenched disorder) could provide 
additional non-trivial interactions among the fluctuating fields. According to 
the arguments presented in this Section, the partition function which describes 
these interactions is analogous to that of the random energy model (REM) of 
Derrida.|12[ Using the known results[|12|] for the REM, it then follows that if 
the interaction, VL-p between local solutions and perturbative fluctuations ex- 
ceeds a critical strength, replica symmetry will be broken. Thus, we identify a 
scenario for obtaining spontaneous RSB from a random spin model. However, 
it remains an open question whether, when the bond randomness is arbitrarily 
weak, the interaction Vl-p exceeds the critical strength necessary for RSB. If 
so, the analogous REM has the 1-step RSB structure and the value of the co- 
ordinate of the step xq is defined by the concrete statistical characteristics of 
the disorder involved. 

The remaining problems as well as future perspectives are discussed in the 
Conclusions. 
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2 THE RENORMALIZATION GROUP AND REPLICA 
SYMMETRY BREAKING 



We consider the p-component </> 4 theory with quenched random effective tem- 
perature fluctuations, which near the transition point can be described by the 
usual Ginzburg-Landau Hamiltonian: 



H[6r,<j>] =Jd D x 



ELi(VMx)) 5 



(2-1) 



+±(r - 5t(x)) ELi tf(x) + \g EL=i 



where the quenched potential 5r(x) is distributed according to Eq. ( |f . 2| ) . In 
terms of the standard replica approach one has to calculate the following replica 
partition function: 



Z n = (JD^expi-HiSr,^})- = J D5rP[dr](J Dfatx) exp{-H[5r, <p]}) n 



J D5t(x) J D(j>^(x) exp 



-±fd D x[Sr(x)Y 



Jd^EUE^mix)} 2 + |[r - 5r(x)] E? =1 E: =1 [0?(*)] 2 



+^el =1 e: =1 [^(^)] 2 [^(^)] 2 ] | 

(2.2) 

where the superscript a labels the replicas. (Here and in what follows all ir- 
relevant pre-exponential factors are omitted.) After Gaussian integration over 
St(x) one gets: 



Z n = JDtf(x)exp 



$d D x\\ e? =1 e"=i m (x)?+ i 2 rj:uj: n a=1 m^w 



+|EL=iE^ =1 ^[C(^)] 2 [^(^)] 2 



where 



9ab = 9$ab ~ u . 



(2.3) 
(2.4) 



To study the critical properties of this system we are going to apply the 
standard RG procedure developed for dimensions D = 4 — e, where e < 1. 



7 



Along the lines of the usual rescaling scheme (see e.g. 0) one gets the following 
RG equations for the interaction parameters g^: 



= £9ab - ^(^9 2 ab + ^9aa9ab + V J2 9ac9cb) , (2.5) 

where £ is the standard rescaling parameter. 

If one takes the matrix g^ to be replica symmetric, as in the starting form 
of Eq. (|2.4j ), then one would recover the usual RG equations for the parameters 
g and u, and eventually one would obtain the well known results for the fixed 
points and the critical exponents 0, |3|. Here, however, we are going to study 
the stability of the RG flows with respect to RSB in the matrix g^. We leave 
the question as to how perturbations out of the RS subspace could arise until 
the next Section, and formally consider the RG eqs. (|2.5| ) assuming that the 
matrix g a b has a Parisi RSB structure (which also include the RS structure as 
a special case). 

According to the standard technique of the Parisi RSB algebra (see e.g. 
P|, |P-0|| ), in the limit n — > the matrix is parametrized in terms of its 
diagonal element g and the off-diagonal function g(x) defined in the interval 
< x < 1: 

9ab -> (g; g(x)) (2.6) 

The RS situation corresponds to the case g(x) = const independent of x. All 
the operations with the matrices in this algebra can be performed according to 
the following simple rules: 



and 



(2-7) 



(9 2 )ab = 9ac9cb -> (c; c(x)) (2.8) 

c=l 



where 



c = 9 2 - Jodxg 2 (x) 

c(x) = 2(g - Jo 1 dyg(y))g(x) - % dy[g{x) - g(y)} 2 



(2.9) 



Using the above rules, and redefining (just to simplify formulae): g a b 
S7i 2 eg a b and £ — ► from the eqs.( |2.5| ) one gets: 



±g(x; = g(x; f) - 4g 2 (x; f) - (4 + 2p)g(x; 0^(0+ 

+2pg(x; Jo 1 dyg{y- £)+p J X dy[g(x; - ^(y; £)] a (2 ' 10) 

^0(0 = fir(0 - ^ + p)~9 2 (0+pl 1 dyg 2 (y,0 (2.11) 

The usual RS equations are recovered if one takes g(x;£) = g>(^) independent 
of x. 

Then, for the fixed point, -^g(x) = 0, = one obtains the following 
equations: 

g(x) - Ag 2 {x) - (4 + 2p)g(x)g+ 

(2.12) 

+2pg(x) Jo 1 dyg(y)+pf* dy[g(x) - g(y)} 2 = 

g - (8 + p)^ 2 + p / c/2/^7 2 (2/) = (2.13) 

JO 



It is interesting to note that the structure of Eq. Q2.12 ) is similar to that for 



the Parisi function q(x) near T c in the SK model of spin-glasses (SG). Although 
the cubic terms in the "order parameter" function g(x) are not included here, 
they could be calculated in the next order of the loop expansion of the RG. 
The most essential difference from the situation in the SK spin-glasses is the 
presence of the term g 2 (x) in Eq. ( |2.12|) (there is no such term in the SK 
model). This is the typical term which is known to produce lstep RSB in 
other SG-like systems ||11|| , and in this case the higher order terms produce no 
qualitative change in the results. 

From the Eqs. (|2.12|) and ( |2.13j ) one can easily find out what should be the 



structure of the function g(x) at the fixed point. Taking the derivative over x 
twice, one gets, from Eq. ( |2.12j ): g'(x) = 0. This means that either the function 



g(x) is constant (which is the RS situation), or it has the step-like structure. 

Let us consider the simplest ansatz, which is the lstep RSB (actually, a bit 
later we are going to argue, that there are no many-step fixed points). Thus 
we assume that 

q(x) = { 9o for < x < x (2 u) 

v ' \ gi for x < x < 1 v ' 

where < Xq < 1 is the coordinate of the step. 

Then, from Eqs. fl2.12p and ( |2.13| ) one gets three equations for three param- 
eters g,g and g^. 
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(4 - 2px )g 2 - 2p(l - x )g 1 g + (4 + 2p)gg = g 



-px g 2 + (4 - 2p + px )g 2 l + (4 + 2p)g 9l = g x (2.15) 

-pxo9o - p(l - xo)^? + (8 + p)<? 2 = § 

Note again, that for g = g\ one obtains the usual RS fixed point equations, 
and the parameter xq drops out from the equations. 

Equations ( 2.15|) have several solutions. Among them there are two fixed 



points which are the usual RS ones 



g = gi=0; 9 = ^ (2.16) 



and 



90 = 5l = -i6FT) ; 9 = W^ { ] 

The first of these corresponds to the pure system, and is stable when the 
disorder is irrelevant according to the Harris criterion. The disorder-induced 
fixed point ( [2.171) is usually considered to be the one which describes the new 
universal critical behaviour in systems with impurities. This fixed point has 
been shown to be stable (with respect to the RS deviations!) for p < 4, which 
is consistent with the Harris criterion since the specific heat critical exponent 
associated with this fixed point is negative and that associated with the pure 
system fixed point is positive. (For p = 1 this fixed point involves an expan- 
sion in powers of (e) 1 / 2 |JJ. This structure is only revealed within a two- loop 
approximation and was therefore not located in the early work 0.) 

However, besides these two RS fixed points there exist the following two 
non-trivial lstep solutions with g ^ g\. 

gu _ 4 ~ P x g 1= g= PX ° (2 18) 

16(px — 1)' 16(px — 1) 



and 



#0 = 0; g x = - 

p(l-30) 



16(p-l)-px (8+p) 

(2.19) 



y 



16(p-l)-px (8+p) 

The stability analysis for the above fixed points with respect to general 
RSB deviations would appear to be technically a rather difficult problem, but 
one can easily check the stability within the considered lstep RSB subspace. 
The calculations, although a bit cumbersome, are straightforward. One has 
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to linearize the "dynamical" equations (|2.10| ) near a given fixed point, and 



then one has to find the eigenvalues. The fixed point is stable only if all three 
eigenvalues are negative. We omit this purely technical analysis and report the 
results only. The pure system fixed point ( |2.16| ) appears to stable (in accordance 
with the Harris criterion) for p > 4, when the pure system specific heat critical 
exponent is negative. 

The traditional RS fixed point ( |2.17] ) appears to be always unstable. The 
three eigenvalues of the linearized equations are: Ai = —1/2, A2 = — § (p~ P f) 

and A 3 = + s (p^i ) • ^ ne instability appears just due to the possibility of creating 
a "step". 



The first lstep RSB fixed point ( |2.18|) appears to be stable in the region of 
parameters: 1 < pxo < 4. 

Finally, the second lstep RSB fixed point (|2.19|) is stable for: 

1 < p < 4; 

(2.20) 

< x < x c (p) = igg^ 

In particular, x c {p = 2) = 4/5 and x c (p = 3) = 32/33 (while x c (p = 4) = 1). 

Note, that in addition to the fixed points listed above there exist several 
other lstep RSB solutions of the eqs. fl2.15D which we do not reproduce here 



because they are always unstable. 

Actually the "physical" fixed point is that in Eq. ( |2.19| ) (with the stability 
conditions ( |2.20| )), and not that in the Eq. fl2.18| ), which has to be considered 
as " unphysical" . The point is that according to the arguments which are due 
to be presented in the next Section, whenever the RSB perturbation comes 
into play, it always requires (according to its physical origin) that g\ < g . On 
the other hand, simple numerical solution of the evolution Eqs. ( |2.10| ) clearly 
demonstrates that if the initial values of the parameters are bounded such that 
gi < go, then, whatever the actual starting values of the parameters are, one 
always ends up at the fixed point ( |2.19| ) and not the one ( |2.18j ). In the opposite 



case, gi > go, the RG trajectories always go to the fixed point ( |2.18 ), but this 
situation must be considered as "unphysical". 

As for the possibility to have many-step solutions of the fixed point Eqs. 
( 2.12 ) and ( 2.13|) , the generalization of the above- considered approach is straight- 



forward, although technically it is much more cumbersome. On the other hand, 
the direct numerical observations of the RG trajectories for the case of several 
steps in the initial function g(x) is still rather simple. Such an analysis convinc- 
ingly demonstrates that whatever the starting conditions are, one always ends 
up in the lstep fixed point (|2.19|) , where the value of the parameter x is the 



coordinate of the "right-most" step in the initial configuration of the function 
g(x). 
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It should be stressed here that, within the present pure RG considerations, 
the value of the parameter xq which defines the fixed points, remains arbitrary. 
In this sense, one can say that in the problem under consideration there exists 
a whole line of fixed points instead of the unique universal one. Formally, the 
value of the parameter xq which characterizes a given fixed point is defined by 
the initial conditions of the RG equations, and in this sense one could argue 
that the critical behavior in systems with such disorder is not universal. 

Another painful problem is to elucidate what is going on if the value of 
the parameter x happens to be in the "instability region" x c (p) < Xq < 1. 
Formally, in this case the RG trajectories go to infinity, and it means that 
one has to take into account the next order RG terms, which hopefully could 
stabilize the situation, just as they do for the Ising (p = 1) case.0 We leave 
these problems for future analysis. 

Consider now what the consequences of the existence of the above lstep 
RSB fixed point are for the critical exponents. The RG equations for the mass 
term of the Hamiltonian in the Eq. (|2.3|) are: 



d n 

—T ab = 2r ab - e[2r ab g ab + p5 ab ^ T cc g cb \ (2.21) 
°^ c=l 

This is a general (one- loop) RG equation for an arbitrary mass matrix r ab . 
In our case the initial matrix is diagonal: T ab (C, = 0) = To5 ab , and it remains 
diagonal (in the higher orders of the RG as well) whatever RSB takes place in 
the interaction matrix g ab . This is just a general consequence of the absence of 
fields which would break the symmetry <p — > —<p. At the lstep RSB fixed point 
(|2.19|) for the rescaling of the (diagonal) mass term one gets: 



2 - e[(2 + p)g - p(l - x )g 1 - px g } 



dim 

(2.22) 



6p(l-xp) 



"l6(p-l)-px (p+8) 

According to the standard scaling relations for the critical exponent of the 
correlation length one finds: 

"M = ^ = ~ + -e , W 1 -* ] (2.23) 

v ; A r 2 2 160- !) -px {p + 8) v ' 

Correspondingly, for the critical exponent of the specific heat: a = 2 — (4 — e)u, 
one obtains: 

/ \ = 1 (4-p)(4-px ) 
a{Xo> 2 C 16(j9- 1) -px (p + 8) { ' ' 
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Thus, depending on the value of the parameter xq one finds a whole spectrum 
of the critical exponents. In particular, the possible values of the specific heat 
critical exponent appear to be in the following band: 

(4 - v) 

- oo < a{x ) < -e± ^ (2.25) 

8(p- 1) 

The upper bound for a(x ) is achieved in the RS limit x — > 0, and it coincides 
with the usual RS result @. On the other hand, as x tends to the "border of 
stability" x c (p) of the lstep RSB fixed point, formally the specific heat critical 
exponent tends to minus infinity. 

Note that, as usual, to obtain the leading fluctuation correction to the critical 
exponent of the correlation functions (usually called rj) one has to study the 
RG fixed points in the next order (~ e 2 ) approximation. 



3 POSSIBLE SCENARIO FOR SPONTANEOUS RSB 

In this Section we will present qualitative arguments showing how RSB per- 
turbations could be spontaneously generated in the random bond model. For 
simplicity we set p — 1. We start by considering the partition function for a 
fixed configuration of Sr(x): 

Z[5t] = J D(j)(x)exp{-H[(f);5T}} , (3.1) 

where 



H\fr 5r] = J d D x ^-[V(j>{x)f + l -[r - 5r(x)]<P(x) 2 + -^{xf) • (3-2) 
for a given realization St(x), the saddle-point equation, 

- A(f>(x) + (t-5t {x)(f)(x) + gct){xf = (3.3) 

has many local minima solutions. We denote such a local solution by ip^\x) 
with i = 1, 2, . . . iVo- If the size Lq of an "island" where 5t(x) > is not too 



small, then the value of ifj^(x) in this "island" should be ~ ±ydr(x)/g, where 
5t(x) should now be interpreted as the value of St averaged over the region of 
size Lq. Such "islands" occur at a certain finite density per unit volume. Thus 
the number of such local solutions, Nq, is macroscopic: iV = kV, where V is 
the volume of the system and n is a constant. An approximate global extremal 
solution §(x) is constructed as the union of all these local solutions without 
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regard for interactions between "islands." Each local solution can occur with 
either sign, since we are dealing with the disordered phase: 



kV 



& a \x ] 5T{x)]=Y J °^ (i) {x) , 



(3.4) 



i=l 



where each Oi = ±1. Accordingly, the total number of global solutions must be 
2 kV . We denote these solutions by [x; St(x)], where a = 1,2,..., If = 2 kV . 
(In this type of symbol we later write simply St for 5t(x).) As we mentioned, 
it seems unlikely that an integration over fluctuations around <p(x) = will 
include the contributions from the configurations of <f>(x) which are near a $(x), 
since <&(#) is "beyond a barrier," so to speak. Therefore, it seems appropriate to 
include separately the contributions from small fluctuations about each of the 
many <&( a )[x; St]. Thus we have to sum over the K global minimum solutions 
(non-perturbative degrees of freedom) $ <Q ) [x; St] and also to integrate over 
"smooth" fluctuations (p(x) around them: 



Z[6r] = fD<p(x) exp \-H 



<S>^[x; St] + ip(x); St 



(3.5) 



/ Dip(x) exp -H[<p; St] x Z\ V ; St] , 



where 



x exp ( - / d D x 



[ X] ST] 2 (p(x) 2 + gQW [x, St]^{x) 



(3.6) 



Next we carry out the appropriate average over quenched disorder. To do 
this, we need to average the nth (n — > 0) power of the partition function. This 
is accomplished by introducing the replicated partition function, Z ni as 



Z n = J DSt J D Va exp (--^/ d D x[ST(x)] 2 - J2 H[<p a ; St]^ x Z n [ Va ;ST] , 

(3.7) 

where the subscript a is a replica index and 
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exp[-EaH[^ aa) }-Jd D xj: a 



(3.8) 

where stands for $[x,<5r]. 

It is obviously hopeless to try to make a systematic evaluation of this repli- 
cated partition function. The global solutions are complicated implicit 
functions of 5t(x). These quantities have fluctuations of two different types. 
In the first instance, they depend on the stochastic variables 5t(x). But even 
when the St(x) are completely fixed, $( Q )(x) will depend on a (which labels 
the possible ways of constructing the global minimum out of the choices for 
the signs {a} of the local minima). A crude way of treating this situation is to 
regard the local solutions ip( l \x) as if they were random variables, even though 
Sr(x) has been specified. This randomness, which one can see is not all that 
different from that which exists in a spin glass, is the crucial one which we 
claim may lead to RSB. Accordingly, we no longer bother to keep track explic- 
itly of the fluctuations in & a \x) due to its dependence on 5t(x). Instead, we 
introduce a distribution function for ip(x). The simplest and the most natural 
distribution function is the Gaussian one: 



P[{iP(x)}] = P[&")[x;5t]] = n.exp -A / d D x[^\x) 2 - ^(x) 2 } 



21 2 



exp(^-± f d D x[& a \x) 2 - <S> {x) 2 } 



(3.9) 



where ipo{x) — \/^ T (x)/g, A is a parameter, which in principle should be defined 
through the original parameters g and u, and <&o(x) is the value of <3> given by 
Eq. ( |3.4| ) when ip^>(x) is replaced by ip®. The above distribution function 
exhibits two Gaussian maxima (with a mean square width equal to A) around 



the values ±v/ 5r{x)/ g in the "islands", over which, on average, 8t(x) > 0. The 

width A reflects the distribution in the magnitude of ip^(x) which results from 
distant "islands" fluctuating between their "up" and "down" states. 

Since the distribution ( |3.9|) is symmetric with respect to the signs of the 
^>W, the term & aa (x)(p a (x) 3 in the partition function (|3.8|) can produce only 



interactions of the order <p , which are irrelevant for the critical properties. 
Therefore, this term can be safely omitted. Using the equation (|3.3|) for the 
energy in a given minimum one easily gets: 
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Then, for the partition function Z n one obtains: 



(3.10) 



J (x) exp ^-^ / (x) 2 - $o(x) 2 ] 2 ^) 



x exp Ug J d D x E a ® aa (x) 4 - \g J d D x E a $ K) WV«W 2 ) • 

(3.11) 

The non-trivial RSB effects we are looking for come from the integration in the 
vicinity of the points ipi(x) = ±ipo(x). Assuming that the parameter A is small 
enough, one can redefine $ a (x) 2 = $o(^) 2 + z a (x). 

Then, for the part of the partition function Z n which contains the integration 
over z a (x) one obtains: 



Z ~ 



IX 



J Dz a (x)exp 




exp 



\9 I d°X Ea Z 2 aa (X) -\ 9 ! d D X Ea *a a (x) [?, Va {xf - $ (x) 



(3.12) 

At this stage the problem (|3.12| ) seems to be similar to that of the Random 
Energy Model (REM) [Jl2| . Leaving apart the details of the rigorous consid- 
eration, one can obtain the correct solution for this problem in the framework 
of the following simplified heuristic procedure. Each term in the exponent of 
the partition function ( |3.12| ) is the sum of n values {^ aa }'s chosen out of the 
total K = 2 kV ones. From the solution of the REM it is seen that the leading 
contribution to the partition function comes from the configurations in which 
in the summations Ea=i z a a one takes n/xo different a a 's, repeated xq times. 
Here xq is a parameter originally defined in the interval 1 < xq < n, which (as 
usual in the replica formalism) turns into 1 > xq > in the limit n — > 0. This 
parameter has to be fixed by extremizing the free energy. If such an extremum 
is achieved for x = 1, then one gets the RS solution, otherwise the system 
appears to be in the lstep RSB state. 

According to the above ansatz for the partition function ( |3.12|) one gets: 



1G 



Z„~(2« y )=olG£ 



n/x 



J Dz ac (x) exp (- ^7 J rf D x4 c (x) 



Here (2 Kl/ ) a: o is the combinatoric entropy factor, and ^7 = ^ — 
After simple algebra one obtains: 



(3.13) 



Z n ~ exp §g 2 A' / dPx E n c L? Et%=i <& 



(3.14) 



- \g* A'x J d D x% E n a=1 + \g 2 A'nxo / c^x$ 2 + f^V In 2) . 
Coming back to the initial problem of integration over the fluctuations (p(x) 



Eq. (|3.7|) , one finds that the second term in the exponent Q3.14Q gives an 
irrelevant shift of the mass term in the Hamiltonian H[ip(x)], while the first 
term in ( 3.14j) provides the RSB structure of the matrix g a b in the interaction 
term: 



1 n 

J 9ab<pl<Pb 
4 a,b=l 



(3.15) 



The matrix g ab appears to have lstep RSB block structure described by the 
parameters (in notations of the Section 2): 



g = g-u- §# 2 A' 



9i 



-u - f# 2 A' 



(3.16) 



9o 



—u . 



Since g, u, and A' are all positive by definition, one finds the following 
restrictions on the values of the interaction parameters: g < 0, g\ < 0, g\ < g 
and ~g > g x . 



In the original REM problem [|12]| , after the partition function is calculated, 
the parameter xq is fixed by extremizing the resultant free energy with respect 
to xq- Here the parameter xq enters into the further problem of integration 
over fluctuations which can be done in terms of the RG procedure. However, 
in terms of the RG technique one usually gets the results only for the singular 
part of the free energy, which is actually small (in r) compared to the whole 
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free energy. All that makes the problem of fixing the parameter xq rather 
non-trivial. 

It should be stressed, however, that aim of the considerations of this Section 
was to identify the physical mechanism which may, in principle, give rise to 
RSB perturbations whose treatment then follows from the usual RG calcula- 
tions of the critical behaviour. It will be very difficult to make our arguments 
really precise. It is clear, in view of the somewhat slippery assumptions made 
above, that the actual contribution to the interactions of the fluctuating fields, 
coming from the non-perturbative degrees of freedom, could have an even more 
sophisticated RSB structure, than the simple lstep form obtained above. How- 
ever, the main point of the present discussion was to demonstrate that, as far 
as the effective interactions of the fluctuating fields are concerned, the RSB 
perturbations could exist in the critical region. To what extend such RSB per- 
turbations are relevant for the critical properties, can then be analysed (as we 
have done) in terms of the traditional RG approach. The results of Section 2 
clearly demonstrate that whenever the weak disorder is relevant (i. e. if the 
specific heat of the pure system is positive), and if spontaneous RSB occurs, 
then the the critical behavior is modified in a dramatic way. 



4 CONCLUSIONS 

We may summarize our conclusions concerning the random bond p-component 
Heisenberg ferromagnet as follows. 

1. The traditional fixed pointsffj, [J in the weakly random ferromagnet for 
the case when randomness is relevant are only stable within the space 
of replica symmetric potentials. Therefore, the corresponding results 
for the critical exponents and other critical properties, are correct only 
as long a replica symmetry breaking does NOT occur spontaneously. 
When randomness is not relevant, our analysis reduces to the standard 
one.§ 

2. Spontaneous replica symmetry breaking has a dramatic effect on the 
renormalization group flows and therefore on the critical properties. 
At first order in e and for p not close to 1 (the Ising limit), the stable 
fixed point corresponds to one step replica symmetry breaking, in close 
analogy with the random energy model. [|l2j Presumably a calculation 



at higher order in e would lead to predictions for how the correlation 
functions would reflect a breaking of replica symmetry. 

At first order in e, there is an instability region near p = 1 where we 
find no stable fixed point in the presence of replica symmetry breaking. 
This result is in close analogy to the known result in the absence of 
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replica symmetry breaking, namely that only at two-loop order does 
one recover a stable fixed point not present in first order. 

The replica symmetry broken fixed point is characterized by a param- 
eter xq whose value is not fixed in the order to which we work here. 
Accordingly, two main possibilities exist. In the first case, the value of 
this parameter may be determined in higher order in e, in which case 
one would have universal exponents, amplitude ratios, etc., as usual 
for a critical point. In the second case, one would have continuously 
variable exponents and a corresponding lack of universality. We have 
no idea what aspects of the randomness, when varied, would lead to 
variation in the critical exponents. 

A key question, which remains unanswered, is whether or not in the 
case of arbitrarily weak randomness our model has spontaneous replica 
symmetry breaking. We have given a scenario by which replica sym- 
metry can be spontaneously broken by interactions, Vlp, between the 
perturbative fluctuations (usually treated within a renormalization 
group context) and fluctuations about local mean field solutions (ig- 
nored in previous treatments). In this scenario we relate the partition 
function due to fluctuations about local mean field solutions to the 



random energy model. |12| The crucial question, which we can not 
answer, is whether Vlp is large enough that the analogous random 
energy model is in its replica symmetry broken phase or not. So, it 
is possible that there is a critical strength in the randomness, below 
which replica symmetry breaking does not occur. 

6. One may mention related work. For instance, the existence of local 
solutions to the mean field equations reminds one of the Griffith phase, 
in which field derivatives of the free energy are anomalous even outside 
the critical region [0|. While the present paper was in preparation 
we learned about similar RSB instability in the RG flows in the 2D 
random field XY model |TJ . 
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